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Abstract 

We consider a medium of multilevel atomic systems interacting with radi- 
ation pulses. A relatively simple technique of analytic calculations is pro- 
posed, which allows revealing all necessary conditions (with sufficient condi- 
tions to be checked separately) imposed on the interaction parameters, for 
which the mean dipole moment of a multilevel atomic medium vanishes, i.e., 
the medium becomes transparent via adiabatic interaction. The proposed 
technique is based on the method of quasienergies and illustrated for three- 
and five-level atomic systems. The necessary conditions for the propaga- 
tion length where the interaction adiabaticity is preserved in the medium are 
obtained. 

Keywords: adiabatic interaction, population transfer, destructive 
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1. Introduction 

Microscopic theory of coherent resonant interaction between radiation 
and individual atoms is sufficiently complete [H, 0, ^ but in case of macro- 
scopic volumes of interaction a number of questions arises which are up today 
open. First, this is the question of spatial scale where the processes occur- 
ring with individual atoms can be realized. Then, this is the question of 
the energy exchange between a pulse and the medium and between different 
pulses travelling in the medium simultaneously. Joint action of nonlinear and 
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dispersive properties of a resonant medium leads to modification of temporal 
and spectral characteristics of an optical pulse and to loss of optical infor- 
mation contained in the incident pulses. Adiabatic processes constitute a 
specific sub-class of coherent interactions, which rely on adiabatic evolution 
of the quantum system 0,0]. As the evolution is slow, the system typically 
remains in a certain laser-controlled quantum state, which is an eigenstate of 
the interaction Hamiltonian. This yields a relatively stable situation. How- 
ever, even at adiabatic interaction the nonlinear dispersion of the medium 
leads in general to parametric broadening of spectra of pulses [6j , which in its 
turn results in rather fast destruction of the interaction adiabaticity • The 
possibility of compensation for the medium dispersion was revealed experi- 
mentally as early as in 1973 in case of medium of three-level V-type atoms [sj]. 
After that the possibility of soliton regime of pulse propagation in such media 
was studied in detail [9j. This mechanism of compensation for the nonlinear 
dispersion of medium is based on destructive interference of dipole moments 
induced by radiation at the adjacent transitions of medium atoms. Another 
striking example of bleaching of resonant medium via adiabatic interaction is 
EIT 0, llll which is based on the so-called dark state formed in a three-level 
system [12]. The dark state does not involve the intermediate state because 
of quantum interference; due to this fact the field-induced dipole moments 
vanish leading to medium bleaching. The purpose of the present work was to 
find the necessary (but not always sufficient) conditions for similar bleach- 
ing of a medium consisting of atomic systems with more complicated level 
diagrams. We describe a relatively simple technique of analytic calculation 
of adiabatic propagation of light pulses and illustrate this technique for a 
five-level atomic system. One of the above-mentioned mechanisms (or their 
combination) underlies the adiabatic bleaching to be considered. We also 
obtain the necessary conditions for the propagation length where the inter- 
action adiabaticity does not break in the medium. 



2. Formalism 

We assume that considered pulses have durations much shorter than all 
times of relaxation and at the same time much longer than the inverse fre- 
quency distance between the closest quasienergies (eigenvalues of the in- 
teraction Hamiltonian), in order to ensure the adiabaticity of interaction. 
Adiabatic propagation of pulses in an atomic medium is described by a self- 
consistent system of truncated equations in running coordinates x and t—x/c; 
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these equations may be represented in the form 
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dE j 



.2tc Nujh dX 



(1) 



dx 



c dE* 



Here A is the eigenvalue of the interaction Hamiltonian, Ej the complex 
amplitude of the pulse field, and other notations are conventional. If the 
system has a quasienergy which remains constant during the overall time 
of interaction, i.e., dX/dE* = 0, then according to flTJ, the induced dipole 
moment in such a system is zero and the pulse propagates in such medium 
without change in shape at the group velocity equal to c. It is known jl5| that 
the quasienergies (eigenvalues) of the system are the roots of characteristic 
equation 



where I is the unit matrix and H the interaction Hamiltonian (in frequency 
units) which may in the resonance approximation be represented in the form 



with o"j j being the projection matrices, f2j the Rabi frequencies of corre- 
sponding transitions i — > i + 1, and the (i-l)-photon detuning (<5o=0) 
(Definition of multi-photon detunings depends on the specific scheme of in- 
teraction, see Appendix). The Rabi frequencies are assumed to be real and 
positive while their phases, which can vary when travelling in the medium, 
are included into one-photon detunings (Aj = Ui + if>i, if > and 

Aj = Uij+i—Ui + cpi, if < 0). By differentiating the characteristic equa- 
tion with respect to E* and setting to zero the derivatives of quasienergy, 
we obtain a system of algebraic equations, from which the needed values 
of interaction parameters may be determined. The problem can, however, 
be simplified. Really, at turning off the fields the interaction Hamiltonian 
is diagonalized and the roots of equation (J2]) go to the following constant 
values: 



The problem of determination of necessary (but not always sufficient) 
conditions for medium transparency is thus reduced to the search of such 
parameters of adiabatic interaction for which the system quasienergy re- 
mains equal to one of resonance detunings always during interaction. By 



det(H - XI) = 



(2) 




(3) 




(4) 
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Figure 1: Several five- level configurations and relevant Rabi frequencies. 



substituting these values successively into equation (j2J) we obtain the needed 
conditions for the interaction parameters. 

For example, in case of a three- level system the root A = is realized 
at exact two-photon resonance 62 = and under the same condition the 
root A = 62 = is realized. The wave function corresponding to these 



quasienergies (dark state) is well known 11 



|^) = cos6»|l) - sin6>|3> (5) 

where tan# = Qi/Q 2 - As the population of the level |2) is zero, the dipole 
moments of 1 — > 2 and 2 — > 3 transitions also vanish. The third root 
A = 5% is realized at degenerate pumping = f2 2 ) and the condition for 
detunings 62 = 25%. As distinct from the previous case, the wave function 
involves all bare atomic states: 

\ifj) = ^sin$|l) + cos$|2) - ^sin$|3) (6) 
v2 ' V2 

where tan$ = \/2fii/<5i.The medium becomes transparent because of de- 
structive interference of dipole moments of transitions 1 — > 2 and 2 — > 3. 
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It should be noted that in case of non-degenerate pumping the root X = 5i 
does not provide transparency, since the dipole moments of the transitions 
1 — > 2 and 2 — > 3 do not equal zero and they do not interfere. So, the 
condition A = 5\ is in this case necessary, but not sufficient. 

Now consider a five-level atomic system and the pulses resonant with 
only the adjacent transitions. Several examples of such diagrams are demon- 
strated in Fig.l. We show in what follows that in such systems six regimes 
of adiabatic bleaching of medium are possible. 

a)Under conditions of exact two-photon resonances ($2=0, #4 = 0) a dark 



state is realized in the system (see, e.g., [l6|) which corresponds to the root 
A=0: 

= — (cos#i cos6> 2 |l) — sin^! cos# 2 |3) + sin^x sin #215)) (7) 
R 

where tan#i = ^1/^2; tan #2 = Vt^/Vt^R = ^/ (1 — sin 2 6*2 cos 2 9i) , This 
superposition does not involve the levels 2 and 4 (destructive interference) 
and the system is well studied in literature for the 



[16|,ll7j (see Fig.2). 

b) The value A = 5i can be realized under conditions 5i = 64 = S 2 /2 
and degenerate pumping at the transition 1 — > 2 — > 3 (Q\ = Q 2 = ^)- 
The wave function does not contain the level 4 and the dipole moments of 
transitions 1 — > 2 and 2 — > 3 are in antiphase leading to their destructive 
interference: 

1*0) = 4( sin <& cos0 2 |l) + V% cos $ cos #2 1 2) - sin$cos6> 2 |3) + sin $ sin 9 2 \ 5)) 
R 

(8) 

where R = \j sin 2 $ sin 2 8 2 + 2 cos 2 6 2 and tan$ = \/2Q,/5i, tan# 2 = ^3/^4. 
This wave function passes to (E} at Qi — 0. For realization of state ([H]) 
the atom should initially be prepared in state |2). Besides, the pulses at 
transitions 3 — > 4 and 4 — > 5 must be turned on in counterintuitive se- 
quence, i.e., 9 2 (t — > —00) — > 0. If also 9 2 (t — > 00) — > 0, the atom 
returns after interaction into state |2).But if Q4 is turned off prior to ^3 
{0 2 {t — > 00) — > k/2), the atom passes into state |5). So, with use of coun- 
terintuitive sequence of pulses at transitions 3 — > 4 and 4 — > 5 with a 
longer pulse at transition 1 — )• 2 — )• 3 (see Fig. 3), it is possible to transfer 
the population from level 2 to level 5. In this population inversion is 

produced between the levels 5 and 4, which may be used for amplification of 
radiation. This state may efficiently be used in atomic level configurations 
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of Fig. lc (W-system). With the same conditions the value A = 84 may be 
realized. 

c) The value A = 81, discussed above, may also be realized under dif- 
ferent conditions, e.g., at 8\ = 82 = 63,64 = 25 1 and degenerate pumping 
of transitions 1 — > 2 — > 3 and 3 — > 4 — > 5 (Q\ = f2 2 and Q3 = f^). 
(Note that the same conditions lead to realization of the values A = 62 and 
A = 63). In this case the wave function involves all bare atomic states, but 
dipole moments of both degenerate transitions are zero (as can easily be seen 
from Eqs (jUJ) below): 

\ip) = — (sin$i sin$ 2 |l) + cos$i sin$ 2 |2) - sin$i sin$ 2 |3)- 
R 

- sin$i cos$ 2 |4) +sin$isin$ 2 |5)) (9) 

where R = a/ sin 2 $1 + sin 2 $ 2 + sin 2 $1 sin 2 $ 2 and tan $1 = tti/5i, tan $ 2 = 
0,3/63. The atom should be initially in state |2) and the pulses turn on in 
counterintuitive sequence. At the end of interaction the atom turns to be in 
state |4). This state may efficiently be used in atomic level configuration in 
Fig. Id for exitation of the hight-lying level. The same state will be obtained 
when Qi = f2 4 and Q3 = f2 2 . 

d) The value A = also may be realized in different ways. For example, 
instead of 62 = 84 = we can require vanishing of three-photon detuning, 
63 = 0, and 62 = 84 7^ 0. In this case, if pumping of the 3 — > 4 — > 5 
transition is degenerate (f^ = f^), we obtain 

|^r) = — (sin $ 3 cos #1 1 1) - sin $ 3 sin0i|3) + cos * 3 sin^i|4) + sin* 3 siii0i|5)) 
R 

(10) 

where R = \J cos 2 $3 + sin 2 Q\ and tan $3 = ^4/^4, tan^i = fii/fi 2 . As 
distinct from the dark state ([7]), it contains four out of five atomic states. 
For realization of this state the atom should be prepared in state |4). 

e,f) Similarly, the value A = 82 can be obtained under conditions 2<5 2 = 
81 + 83 = 84,0,1 = ^4, 1^3 = 2 and the value A = 83 under conditions 
81 = 8 2 , 38 3 = 82 + 84, Qi = Q 2 - 

3. Propagation 

Although in adiabatic approximation all the above- considered schemes 
lead, as follows from Eq.(jTJ, to adiabatic bleaching of the medium, we cannot 
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Figure 2: Dynamics of atomic level populations (pjj) in M-system and the sequence of 
pulses. Shapes of pulses are chosen to be Gaussian and the value of the one-photon 
detuning is A=10/T. 

determine the criterion of interaction adiabaticity for a single atom knowing 
only a single value of quasienergy. Since we are interested in the propaga- 
tion problem, we must determine the adiabaticity criterion for the medium 
and it may be found directly from the propagation equations. By combining 
the Schrodinger equation with truncated Maxwell equations, we obtain for 
a medium consisting of M-type atoms (Fig. la) the following equations de- 
scribing variations of one-photon detunings (frequencies of pulses) and pulse 
intensities (in running coordinates x, r = t — x/c) during propagation in 
medium: 

dn\ d\h\ 2 



dx ^ l dr 

dn 2 2 <9(N 2 + N 2 ) 



-92 



dx '' dr 

M = ^ /(N 2 + N 2 + N 2 ) (n) 

dx dr 

gOj = <9|& 5 | 2 
dx ^ 5 dr 
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dAj _ d Re(b*b i+l ) 
dx 1 dr £li 

where = 27iNui\d iti+1 \ 2 / he (d^ i+ i is the matrix element of the dipole mo- 
ment of respective transition). For other atomic level schemes equations for 
intensities differ in only the signs of derivatives with respect to populations. 
For example, in case of ladder (cascade) configuration we have to change 
the signs in the second and fourth equations for field intensities (TTTT) . The 
choice of sign is determined by whether the field is absorbed or emitted in 
the specific transition. If the pumping is degenerate (i.e., if the pulses at 
adjacent transitions have the same frequency), the right-hand-sides of cor- 
responding equations should be added. Since it is the time derivatives that 
stand in rhs of these equations, by use for atomic populations the expressions 
following from the wave functions obtained above, we take into account the 
first nonadiabatic correction which lead to the change in spectral and tempo- 
ral characteristics of pulses. Correspondingly, the conditions of smallness of 
these changes are criteria of adiabaticity in the medium. We will now demon- 
strate this statement for a traditional STIRAP-chain in a M-system and for 
a W-system (whose dynamics of travelling in the medium is not studied to 
the best of our knowledge). 



4. Traditional STIRAP-STIRAP-chain in M-system. 

We will restrict ourselves to the case of equal strengths of adjacent transi- 
tions and to the conditions Q± = Q3 and ^2 = ^4- With use of wave function 
(J7|) we obtain directly from Eqs. ffTOj) 



dx dr 

do 2 



<Z7r(N -N 2 ) 



dx dr 



fc(N 2 - N 2 ) (12) 



As follows from the last equation, the self-phase modulation in the medium 
is absent, i.e., the conditions for multiphoton detunings given at the medium 
entrance are preserved during propagation. The quantity + Q 2 , = Qq 
is the integral of motion. Hence, the entire dynamics of propagation of 
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pulses in medium is described by the angle 8(x,t). Since f2 2 = Q 2 l sin 2 9 and 
Q 2 , = JIqCos 2 ^, we finally obtain for the mixing angle 9(x,t) the equation 

£ + *0^ = O (13) 
ox IIq t 

where f{9) = (1 + 2 cos 2 # sin 2 9)/(l - cos 2 9 sin 2 9) 2 . As follows from the 
obtained equation, the variation of the angle 9(x,t) can be neglected at the 
lengths of propagation satisfying the condition 

W m K W K<1 (14) 

This is just the necessary criterion of the adiabaticity of interaction in the 
medium, so that at the travelling lengths satisfying the condition f|T4|) the 
influence of the medium on the process of coherent population transfer may 
be disregarded. Equation f|T3|) may be solved analytically also without this 
limitation. The solution determined by, e.g., the method of characteristics 
may be represented as 

9(x,t) = 9 {t-x/u) (15) 

where 9q is a function given at the medium entrance and u the nonlinear 
group velocity defined by the equation 

[ T n 2 dt = qxf(9) (16) 

J T — x/u 

Because of the nonlinear group velocity the time derivative of the function 
9 will increase infinitely at the travelling lengths qxf(9)/Vt^T ~ 1 which 
leads at this lengths to breaking of interaction adiabaticity. Note that in the 
regime linear with respect to probe field (small mixing angles) the function 
f(9) goes to unity and the shock wave does not form which allows using this 
system for information storage. 

5. Population transfer in W-system 

For a W-system in case where the wave function (JTj) is realized, we obtain 
for intensities and frequencies of pulses the following system of equations 
(transition strengths are assumed, as above, to be equal): 

dQ 2 _ 2 
dx dr 
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Time (units of T) 



Figure 3: Dynamics of atomic level populations (fijj) in W-system and the sequence of 
pulses. Shapes of pulses are chosen to be Gaussian and the value of the one-photon 
detuning is A=10/T. 

«S = _«S= - N . (17) 

ox ox dr 

OA, 



dx 







The self-phase modulation is, as in previous case, absent, since in transition 
2 — > 1, 3 (V-system) conditions of compensation for dispersion are met and 
in transition 3 — > 4 — > 5 (inverted A -system) the conditions for EIT are 
satisfied. The quantity 0| + f2 2 = Q% is the integral of motion. Using the 
definition of mixing angles 82, we obtain directly from ( IT7|) the following 
equation: 

~dx~ + ^^r - QlA ^r~ (18) 

where u± is the nonlinear group velocity, always lower than c, which is defined 
as 

1 1 2g^ 2 (2 + tan 2 g 2 ) 

Ul c (5 2 + fif(2 + tan 2 2 )) 2 1 ' 
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The rhs of Eq.f llSj) is caused by depletion of the field Q\ via transformation 
of photons at Qi into photons of the field f^. The factor A in f llSp equals 



A = 2q51 (20) 

(S* + nl(2 + tan 2 e 2 )Y 1 UJ 

We see that at the lengths in the medium satisfying the inequality 

qxT/(5iT) 2 « 1 (21) 

either of these effects is negligibly weak. Before turning on the field Sl 3 (i.e., 
at 02 = 0, see Fig. 3) the solution of the Eq. (|T7|) is 

ttl(t,x) = nl^t-x/ux) (22) 

Inequality ( 12"T|) has in this simple physical meaning: the group delay 

of the pulse in the medium is much shorter than its duration. After turning 
off the field ^4 (i.e., at 82 = vr/2) we again obtain the same solution. At the 
lengths 4qxQl/5fTi ~ 1 the nonlinear group velocity will lead to formation 
of the shock-wave front, i.e., to breaking of interaction adiabaticity {T\ is 
duration of Q4 pulse). The full dynamics of propagation of the field ft± requires 
joint solution of the equation ( fl8i) and the equation for the mixing angle 62'- 

902 Aq f(f) ^_qBdS^ 

^-^MM^--^^r (23) 



where 



/l(02,*) 



cos 80 sin 2 $ 



(2cos 2 8 2 + sin 2 8 2 sin 2 $) 2 



B = stn28 2 cos 2 $ 1 m) 

(2cos 2 ^ 2 + sin 2 ^ 2 sin 2 $) 2 ((5 2 + 2fi 2 ) 1 ] 

The travelling length where the variation of the angle 82 may be neglected 
should satisfy the condition 

qxTJinoT,) 2 « 1 (25) 

Denominators of expressions (|2~T!) and (25) provide criteria of interaction 
adiabaticity for a single atom, while the whole expressions give criteria of 
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adiabaticity for the medium, i.e., they determine the travelling lengths where 
the propagation effects may be neglected. The medium parameter qxt may 
be expressed in terms of center of line absorption coefficient ao = q/T, qxT\ = 
ocqxTTi. Then, under condition r << 8\ the expression (|2ip may be rewritten 

as 

otoT 2 x 
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« TT « 1 



which means that the linear absorption in medium, caused by the atomic 
linewidth, is negligibly small. However, by integrating the first equations 
([IE]) over time we obtain for the pulse energy Wi = f°° Q 2 dt : 

dWi ni l2 , , x „ |2 



dx 



qx(\b 2 \ 2 (+oo) - \b 2 \'(-oc)) 



In case of complete transfer of population we have |&2| 2 (+oo) — |&2 1 2 ( — oo) = 
— 1. Hence, Wi(x) = —qix + Wo, with Wo being the pulse energy at the 
medium entrance. So, the length at which q±x = Wo is the length of full 
depletion (absorption) of pumping radiation . So, at the propagation length 
satisfying the condition 

qx qx 

< -£=r « 1 



we obtain full absorption of the field Oi caused by the presence of the field 

The main distinction of the equation ( 123]) describing the dynamics of 
propagation of the mixing angle 9 2 from the similar equation ( 1T5]) for the 
M-system, is the fact that the nonlinear velocity of propagation of angle 
6*2 exceeds the speed of light, so the population transfer in this system is 



superluminal [18(|. This situation is similar to that of superluminal population 



transfer in a tree-level system via bright state considered in [18] . Really, by 
neglecting the rhs of Eq.( ]2"2"]) in case of large one-photon detuning, we obtain a 
solution similar to that obtained above for the M-system, but with replacing 
x by -x: 

9{x,t) = e (r + x/v), I tt Q dt = -2qxf l {d 2 ^) (26) 

J T+x/v 

More detailed analysis of propagation of pulses with taking into account 
the first nonadiabatic correction requires combined solution of the system of 
equations ffTB"]) and f l2"B"]) ; it is out of scope of the present work. 
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7. Conclusion 

We have developed and analyzed, for the example of a five-level system, 
a relatively simple method for finding necessary conditions under which adi- 
abatic short light pulses can travel in the medium without distortion of the 
shape and phase. It was shown that the necessary conditions for such bleach- 
ing of medium (i.e., vanishing of dipole moments induced in the medium at 
the frequencies of all interacting fields) is the equality of system quasiener- 
gies to one of the values of single- or multi-photon detunings of resonances. 
In this case the medium bleaching may be caused by both interference of 
quantum states and interference of dipole moments of different transitions 
if pumping is degenerate. The existence of such regimes means that at the 
propagation lengths where the interaction adiabaticity does not break, stable 
superposition states are produced all over the medium and these states may 
easily be controlled by adjusting the parameters of pulses. Finding of the 
sufficient conditions requires determination of roots of algebraic equations 
which may be done only numerically for the order of equation higher than 
three. However, the sufficiency of found necessary conditions may easily be 
checked directly by calculation of the induced dipole moments. In addition, 
we demonstrated for the examples of M- and W-systems the possibility to 
find a generalized criterion of adiabaticity of interaction in the medium with- 
out determination of all eigenvalues of interaction Hamiltonian; we analyzed 
population transfer in these systems and its peculiarities during propagation 
in the medium. 
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